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ABSTRACT 

The Sunyaev-Zeldovich (SZ) effect gives a measure of the thermal energy and electron pressure in groups 
and clusters of galaxies. In the near future SZ surveys will map hundreds of systems, shedding light on the 
pressure distribution in the systems. The thermal energy is related to the total mass of a system of galaxies, 
but it is only a projection that is observed through the SZ effect. A model for the 3D distribution of pressure 
is needed to link the SZ signal to the total mass of the system. In this work we construct an empirical model 
for the 2D and 3D SZ profile, and compare it to a set of realistic high resolution SPH simulations of galaxy 
clusters and groups, and to a stacked SZ profile for massive clusters derived from WMAP data. Furthermore, 
we combine observed temperature profiles with dark matter potentials to yield an additional constraint, under 
the assumption of hydrostatic equilibrium. We find a very tight correlation between the characteristic scale in 
the model, the integrated SZ signal, and the total mass in the systems with a scatter of only 4%. The model only 
contains two free parameters, making it readily applicable even to low resolution SZ observations of galaxy 
clusters. A fitting routine for the model that can be applied to observed or simulated data can be found at 
|http : / /www . phys . au . dk/~ haugbo el/ software . shtml} 
Subject headings: 



1. INTRODUCTION 

When Cosmic Microwave Background (CMB) photons 
pass through galaxy clusters, they Compton up-scatter on hot 
electrons, making a small increment (decrement) above (be- 
low) the peak of the CMB primary spectrum. The size of 
this distortio n in the CMB spectrum, the Sunyaev-Zeldovich 
(SZ) effect dSunvaev & Zeldovichlll972T) . is proportional to 
the electron pressure integrated along the line-of-sight. Be- 
sides being an important probe of the physics of the intr- 
acluster medium, the SZ effect is a promising tool in cos- 
mology: The signal produced by a cluster is practically red- 
shift independent and can be observed t o high redshifts. The 
Planck surveyor satellite (Tauber 2000) will produce a clus- 
ter catalogue with up to ^10.000 (suffi ciently massive) clus- 
ters out to z ~ 1 (Schaf er et al.l l200"6h . By combining SZ 
and X-ray observations of relaxed clusters the Hubble con- 
stant Hq can be derived, limits can be put o n f2 m , and the gas 
mass fraction in clusters can be measured dBonamente et al.l 
l200l lLaRoque et all 12006b . The SZ effect was first de- 
tected in three ne a rby cl usters more than two decades ago by 
iBirkinshaw et al.l d!984l) . but the last five years, due to ad- 
vances in sub millimetre receiver technology, routine mea- 
surements have been done for 0(100) clusters using facili- 
ties such as the OVRO and BIMA telescopes. Right now a 
large number of telescopes have started observing the SZ sig- 
nal or are under construction (e.g. ACBAR, CARMA, SUZIE 
III, SPT, APEX etc.), and in the near future ALMA will 
gradually become online enabling unprecedented resolution 
and sensitivity for making detailed observations of individ- 
ual clusters. Current SZ surveys of galaxy clusters (see e.g. 
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iBonamente et al.ll2006l for a set of current observations) have 
observed the unresolved integrated SZ signal from clusters, 
but with CARMA, ALMA and the SPT also the SZ signal as 
a function of radius fr om the cluster centre will be measured. 
lAfshordi etaf] d2007l) has used the 3rd year WMAP data to 
extract SZ images from 193 clusters. Stacking them they have 
obtained an averaged SZ profile for clusters with Tx > 3keV. 

Spatially resolved SZ images probe, in a way complemen- 
tary to X-ray observations, the distribution of thermal energy 
in clusters which in turn is closely linked to the total cluster 
mass. Furthermore, since the SZ effect essentially depends 
on the electron pressure, the physics going into understanding 
it is simpler and more robust, than is the case for the X-ray 
emission. 

In this paper we use realistic high resolution N-body/SPH 
simulations of galaxy clusters and groups together with ob- 
served temperature profiles of nearby groups and clusters to 
predict the radial SZ profile of different types of systems of 
galaxies, and compar e our results to the averaged profile of 
lAfshordi et al.l (120071) . We introduce a universal fitting for- 
mula, with only two free parameters, that can be employed 
in future observations. Furthermore, we show that fitting this 
profile enables a precise estimate of the total mass in the sys- 
tem. In the next section we describe the computer experi- 
ments, that are used to construct the synthetic SZ profiles. In 
section 3 we discuss the simulated data and present our re- 
sults, and in section 4 we discuss and provide our conclusions. 

2. COMPUTER EXPERIMENTS 

We use 24 galaxy group and galaxy cluster simulations to 
study the SZ effect for systems of virial temperatures of about 
1 to 6 keV. The models include 12 groups of approximately 
the same mass (M VI > = 0.8-1.1 x 10 14 M Q ), 1 1 identical clus- 
ters with M vir = 2.7 x 10 14 M Q ("Virgo" clusters), but simu- 
lated with different gas physics, and a single large cluster 
of M vir = 12 x 10 I4 M Q (a "Coma" cluster). The models are 
re-simulated from a low resolution cosmological dark matter 
simulation, where the halos are identified with a halo finder. 
The particles are traced back in time to a initial redshift Zmitiah 
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and the virial volume is repopulated with both gas and dark 
matter at a high resolution. The code includes radiative cool- 
ing, star formation, supernova feedback, chemical evolution 
and back reaction from a redshift dependent UV field. The 
different models are summarised in Table Q] For a detailed 
description of the different gas physics going into the code in 
general , and the "Virgo" clu ster simulations in particular we 
refer to lRomeo et al.l (I2006I) . The groups were all selected at 
random, the only criterion being their virial mass, and there- 
fore they comprise a cosmological fair sample of groups with 
M Vir = 0.8-l.l x 10 14 M o . 

The different masses of the clusters allow us to probe the 
mass dependence of the SZ profile, while the "Virgo" mod- 
els are used to investigate how robust our predictions are for 
the SZ profile with respect to assumptions on the underly- 
ing physics. The groups, being a statistically unbiased sam- 
ple, give limits on cosmic variance for the given virial mass 
(1O 14 M ). 

We have divided the groups into three different classes ac- 
cording to their morphology and evolutionary history: Groups 
with merging activity (186,231,239,262), fossil groups 
(189,228,236,244) and normal groups (190,233,247,276). A 
group is considered fossil if the difference in apparent R-band 
magnitud e of the first and s econd brightest galaxy is greater 
than two dJones et al.ll2000l) . and a group is merging if it by 
visual inspection has significant merging activity in the core, 
or if the rms scatter in 3D radial shells of the temperature, 
pressure and density, is comparable to the average value in 
the shell. 



Name T x [keV] 7;„ g [keV] Af 18o [lO' 4 M ] Comments 



Groups 




Re-simulated groups 


with Zi„i,ial = 39 


186 


0.99 


1.08 


0.92 


Merging group 


189 


1.07 


1.09 


0.91 


Cool-Core, Fossil group 


190 


1.21 


1.22 


0.99 


Normal group 


228 


1.07 


1.10 


0.85 


Cool-Core, Fossil group 


231 


0.99 


1.08 


0.84 


Merging group 


233 


1.03 


1.09 


0.82 


Normal group 


236 


1.00 


0.99 


0.73 


Fossil group 


239 


0.91 


1.11 


0.81 


Merging group 


244 


1.13 


1.08 


0.80 


Fossil group 


247 


1.01 


1.15 


0.90 


Normal group 


262 


0.97 


1.03 


0.88 


Merging group 


276 


1.02 


1.07 


0.82 


Cool-Core group 


Clusters 




Re-simulated clusters with Zinitial = 19 


AY-SW 


2.13 


2.07 


2.32 


Arimoto-Yoshii IMF 


AY-SW-8 


2.18 


2.10 


2.36 


8 times resolution 


AY-Vol39 


2.18 


2.18 


2.36 


Zinitial — 39 


AY-SWx2 


2.18 


2.04 


2.40 


2 times SNII feedback 


AY-SWx4 


2.17 


2.09 


2.38 


4 times SNII feedback 


AY-PH0.75 


2.23 


2.10 


2.41 


preh. 0.75keV/part @ z=3 


AY-PH1.5 


2.27 


2.14 


2.37 


preh. 1.5keV/part @ z=3 


AY-PH50 


2.19 


2.07 


2.41 


preh. 50keV ■ cm 2 part @ z=3 


AY-COND 


2.26 


2.08 


2.43 


Thermal conduction 


Sal-SW 


2.24 


2.10 


2.42 


Salpeter IMF 


Sal-WFB 


2.27 


2.23 


2.41 


Weak feedback 


Coma 


5.57 


5.27 


10.4 


Relaxed massive cluster 



TABLE 1 

Characteristics of simulated clusters 



3. ANALYSIS 

The SZ profiles of the simulated systems can only be used 
as templates for an universal fitting formula, if the profiles 



are in accordance with observations. In Fig. Q] are shown the 
spherically averaged SZ pr ofiles compared to th e only cur- 
rently published SZ profile dAfshordi et al.l d2007l) ). obtained 
using data from the WMAP satellite. The observed points 
may indicate a slightly more bended profile than the simula- 
tions, but within the 1-a error bars there is fairly good agree- 
ment. From the simulations there is a clear trend towards 
steeper profiles, for more massive and relaxed clusters and 
groups reflecting differences in the underlying DM potentials 
(see Fig.|2|. This should be recalled when constructing an av- 
erage observational profile, because the average profile may 
not represent a true ("universal") physical profile, but rather a 
smeared average of the real profiles. 

3.1. The inner part of the Sunyaev-Zeldovich profile 

The central core of the SZ profile (r < 0.25 r^o cannot be 
probed by WMAP, due to its limited resolution of at most 
0.12°. To extend the dynamic range of the observations we 
combine gravitational potentials from the simulations with 
observed average X-ray temperature profiles. Under the as- 
sumption of hydrostatic equilibrium in the core of the sys- 
tems, we can then predict the inner part of the SZ profile. 

The Compton v-parameter, which determines the overall 
temperature decrement in the CMB radiation due to the SZ 
effect, is proportional to the integrated pressure P of the elec- 
trons along the line of sight 

f f <J T P(l) f k B T e 

y = y(l)dl= / - L Ardl= / a T n e -^dl. (1) 

J J m e c L J m e c L 

Assuming hydrostatic equilibrium and spherical symmetry 



we can use the proportionality of P, and the Compton y- 
parameter in a volume element, y(l(r)), to obtain 

ln ^ =ln m = -^J ro w^ dr (3) 

We have constructed three different averaged gravitational po- 
tentials based on the normal groups, the fossil groups and 
the "Virgo" cluste rs. For the temperature profiles we use ob- 
served clusters by [Vikhlinin et al. (2005), and make one av- 
erage profile constructed from the galaxy groups with virial 
temperatures less than 2.5 keV, and one constructed from the 
massive clusters in the sample with virial temperatures be- 
tween 3.5 and 8.5 keV. The dark matter distribution, derived 
from the simulations, is quite robust against the specific gas 
physics involved in the simulation. This is demonstrated in 
Fig. where it is seen that all the "Virgo" models, in spite 
of the different gas physics, have essentially the same gravita- 
tional potentials from 0.02 riso and outwards. However, there 
are systematic differences between different types of systems 
with the same mass. Fossil groups are more relaxed com- 
pared to normal groups, and their mass distribution is similar 
to relaxed clusters, like the "Virgo" and "Coma" models from 
0.02 r^o and outwards (see Fig.|2|i. Therefore the above re- 
construction procedure gives a good idea of future resolved 
SZ observations of the cores of clusters of galaxies, and our 
six resolved profiles are a fair sample of what can be expected. 

3.2. Temperature profiles 

From Fig.Q]it is clear, that albeit the SZ profiles constructed 
from the observed temperature profiles and the dark matter 
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FIG. 1. — Spherically averaged 3D SZ profiles computed from the simulations, the stacked profile of Afshordi et al. (2007), and derived 3D profiles from 
observed X-ray temperatures combined with average gravitational potentials. The profiles are normalised to the critical density Pcrit> antl the typical temperature 
of the system, Tx, defined as the average projected spectral-like temperature in the interval 0.1 r^oo < r < 0.4 r2oo- To reduce the visual scatter only average 
profiles from the simulations are shown. 
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FIG. 2. — Dark matter potentials measured as the equivalent rotational 
speed, both for specific groups and clusters, and for the average "Virgo" clus- 
ter, normal and fossil group. The curves are normalised with the virial mass. 
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FIG. 3. — Mass-weigthed 3D temperature profiles for the simulated systems 
compared to observed an average group profile (1.6 keV< Tx < 2.6 keV), and 
an average cluster profile (3.5 keV < Tx < 8.5 keV). 
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density profiles are in approximate agreement with the pro- 
files extracted from the simulations, they have a tendency to 
be more peaked at the centre. This can be traced to differ- 
ences in the observed and simulated temperature profiles (see 
e.g. Eq. [3), and to the typical masses of the observed sys- 
tems. The different temperature profiles are normalised with 
a global temperature Tx- To enab le direct compar i son to the 
observed temperature profiles by Vikhlinin et al.l d2005l) we 
compute T x as the average projected spectral-like temperature 
dRasia et alj|2005l) in the interval 0.1 r 2 oo < r < 0.4 r 2 oo- All 
profiles converge to the same universal curve for r > 0.2riso, 
but in the inner part of the systems there are differences (see 
Fig-Ell. The simulated systems have all flat profiles at the core, 
which are in good qualitative agreement with observations of 
clusters of similar virial mass, but the average observed pro- 
files have lower normalised temperatures towards the centre, 
with more or less the same offset between the "Coma" clus- 
ter and the observed cluster profile (constructed from clus- 
ters with 3.5 keV < Tx < 8.5 keV), and between the average 
simulated and observed group temperature pr ofiles (the latter 
constr ucted from groups with Tx < 2.5 keV IVikhlinin et all 
d2006l) ). The offset compared to observations and relatively 
flat temperature profiles are well known problems for simula- 
tions invoking radiative cooling and feedback processes (e.g. 
iBorgani et al.l2004t|Pratt et alj2006b . but it is further accentu- 
ated by an offset in mass between the observed and simulated 
groups and the observed clusters, and the simulated "Virgo" 
clusters. 

3.3. A universal SZ profile 

As can be seen in Fig. [TJ all the SZ profiles have nearly the 
same form, but the overall normalisation, tension and slope 
of the different profiles depends on the mass, and the specific 
cluster. To construct a simple model we have tried to fit a 
variety of combinations of beta profiles and exponential pro- 
files for the density combined with a polytropic or isothermal 
equation of state, but it does not yield a satisfactory fit. The 
correct temperature to use, when constructing a SZ profile, is 
the mass weighted, and it does not necessarily agree with the 
tempe r ature inferred fro m X-ray observations dBonaldi et al.l 
d2007l) : lAfshordl d2007l) ). which may explain why the above 
combination works well for X-ray observations, while not so 
for our sample of SZ profiles. 

To circumvent this problem we are using a novel universal 
profile that directly fit the data, without assuming anything 
about the underlying temperature and density distributions, 
while using as few parameters as possible. The main morpho- 
logical differences are in the slope/tension and in the normal- 
isation of each profile, and it is indeed possible to construct 
a model, that with only two free parameters, a normalisation, 
and a characteristic scale can fit the full set of simulations and 
observations in detail. 

The profiles extracted from the simulated data do not per 
se contain any errors, but a measure of the natural scatter in a 
profile y(r) at a given radial distance, is the variance er(r) of 
y(r) inside the radial bin at r. We have used this variance esti- 
mate to construct a goodness of fit parameter \ 2 for each pro- 
jected 2D and spherically averaged 3D profile y2o(f), y3D(r) 
given by 



X2D,3D [ 



1 



y2D,3DW-3S, i3D (r,p) 



'2D.3D 



(r) 



■dr, (4) 
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FIG. 4. — The goodness of fit parameter for the different models. It is 
mostly one group with extreme merging activity, and the very smooth models 
for the "Virgo" cluster, with thermal conduction or entropy floor, that the 
model have problems with fitting. 



nal scatter in relaxed systems is much smaller, giving stronger 
constraints than from merging systems. This is sensible, since 
the observational scatter among relaxed clusters is smaller 
too. 

We have found that using a triple exponential model gives 
an excellent fit to the data: 



jSdO) = y m [/i exp(-r/r )+/2exp(-r/(a 2 r )) 
+/3exp(-r/(a 3 r ))] , 



(5) 



where f +/ 2 = 1. While this model at first sight might appear 
complicated, it has two attractive features: (i) The 2D profile 
can be derived in a closed form from the 3D profile, (ii) Most 
of the parameters can be fixed to global values, leaving only 
two free parameters for fitting all systems considered in this 
paper. With only two parameters in the model it can readily be 
applied to future observations, even if they only map the SZ 
profile with a few observational points. By integrating along 
one of the axes we find the 2D profile to be 



y£,(r) =2Y m [/i^i(r/r )+/ 2 ^ 1 (r/(a 2 r )) 
+f3K 1 (r/(a 3 r ))] r, 



(6) 



where K„(x) is the modified Bessel function of the second 
kind. Minimising x\d simultaneously for all models over the 
range 0.01 r\$o < r < 1 .9 rigo, while varying fi and a,- as global 
parameters, and Y m , ro for each model, treating the x—, y—, and 
z-projections as different clusters, we find the global best fit 
parameters to be 



fx =0.043 
/ 2 = l-/i= 0.957 
/ 3 =5.9 
a.2 = 0.36 
0-3=0.115, 



(7) 
(8) 
(9) 
(10) 
(11) 



where y™ D 3D (r, p) is the model with parameters p. The inter- 



In Fig. |4] we see how these yield very reasonable \ 2 f° r a ll 
projections, and the 3D profiles. The only systems that have 
a x 2 significantly larger than one are either merging, or char- 
acterised by a very smooth profile with almost no variation in 
the radial bin, and hence a very small <r(r) (see e.g. Fig. [5]). 

The two free parameters, the overall scale ro, and the overall 
normalisation Y m scale with r^o and with Y\%q, the integrated 
SZ effect inside riso, respectively, but there is a large scatter 
between different clusters, where the tightest relation is found 
for the fossil groups, that act as an extension of the clusters 
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FIG. 5. — Examples of the model fits to the SZ profiles with error bars. 
Two of the profiles represent some of the poorest fits, and illustrate why they 
are poor fits: The 239 group has a poor fit because of the merging activity, 
the flatness inside- and spike at O.lriso is due to the two central galaxies, 
separated by approximately 0.4>'i8o, while the AY-Cond "Virgo" cluster has 
a poor fit because the thermal conductivity in the model gives a very smooth 
cluster with small error bars. The two others are a representative group and 
cluster, both with well fitted profiles. 





FIG. 6. — The fitted profiles for the observed average profile for massive 
clusters (top), and for the constructed profiles assuming hydrostatic equi- 
librium, and using average dark matter profiles from the simulations, and 
observed temperature profiles for groups (bottom left) and massive clusters 
(bottom right). The red, yellow, and blue curves are for the average "virgo" 
cluster, fossil group, and normal group DM potentials. 

with almost the same scaling relation. They are related as 



=15 



r =280 



'180 



10" 6 Mpc- 
lOOOkpc 



0.9 



10" 6 Mpc 2 



1/2 



kpc 



(12) 



(13) 



To push the boundaries of our fitting formula for other sys- 
tems we have also applied it to the lAfshordi et alj d2007l) 
data, and the profiles for the centre of the systems, derived 




FIG. 7. — The slope r ( ) (top) and the overall normalisation (bottom) in the 
model. The fossil groups act as an extension of the clusters, with the same 
scaling relation. 

from observed temperature profiles (see Fig. |6). The aver- 
age SZ profile is described well by the formula , while there 
are problems with fitting the core (r < 0.05 rigo) of the cluster 
set of central profiles (the lower right panel in Fig. [6]). Even 
though some of these combinations (e.g. a normal group DM 
potential combined with a massive cluster temperature pro- 
file) are extreme, the lack of agreement may be because in 
the central parts of the observed systems there are signifi- 
cant non-the rmal contributions to the pr essure balance from 
e.g. an AGN (|Roychowdhury et al.l2005l) or from cosmic rays 
dPfrommer et al.ll2007h . and hence the assumption of hydro- 
static equilibrium does not apply. Or it may be that the sim- 
ulated systems do not include an adequate description of the 
physical processes in the centre of the systems. Nonetheless, 
the contribution of the central 0.05 r\%o to the total SZ signal is 
approximately 5%, and the parameters are not much affected 
even if we cannot reconstruct the innermost part of the profiles 
with perfection. 

3.4. Estimating the total mass in a system 

The SZ profile of the systems seems to be universally well 
described by only two parameters, except for possibly in the 
central parts of the clusters. It measures the distribution of 
thermal energy in the system, and is therefore related to the 
gravitational potential, if we assume the system is relaxed. 
Using assumptions about hydrostatic equilibrium and using a 
phenomenological approach sever al "Fundamental plane" re- 
lation s have been constructed (e.g. lAfsho rdi 2007; Ver de et aT] 
12002b . The main ingredients for the Fundamental plane has 
been the integrated SZ effect and a characteristic scale in the 
system, for example Rsz.i, the radius enclosing half of the 
integrated signal. This has given a relation between the in- 
tegrated SZ effect, a charact eristic scale, an d the total mass 
with a rough scatter of 14% dAfshordill2007l) . With our fit to 
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FIG. 8. — The estimated mass (from Eg. 114) against the real mass of the 
systems (top), and the relative error (bottom). The dotted line indicate the rms 
scatter. For compar ison we have inc luded the error on the estimated masses 
using the method of Afshordi (2007) as light grey symbols. 



the profiles we get a very precise measurement of this char- 
acteristic scale, that take into account the different bending of 
the profiles. Fitting the total mass as a function of Y^o and r$ 
we find (see Fig.0 4 



M m = 1.02 x 10 14 M G 



'180 



0.57 



ro 



0.27 

(14) 



\Q- b Mpc 2 J \1000kpc 

with only a 4% scatter. In lAfshordil (120071) a phenomeno- 
logical model is used to construct a set of, essentially, one 
dimensional models for the clusters. It includes normal or 
tophat distributed controlling parameters, with broadening in 
agreement with observed clusters. To check if the smaller 
scatter we see for Eq.[14]is due to more similarity in the sim- 
ul ated systems, we have tried to use the fundamental plane 
of lAfshordil (120071) . on the simulated systems. We find good 
agreement with a 15% scatter for our data (see Pig. [5J. 

The first resolved SZ profiles will be obtained for nearby 
massive clusters. To make a crude test for how our fitting for- 
mula works on real data, we have made a set of mock observa- 
tions taking as a starting point the most massive cluster in our 
sample, the "Coma" cluster. We use five logaritmically spaced 
rings at a distance of [0.09,0.18,0.35,0.7, 1.4] r im . The reso- 
lution of the first generation of SZ instruments is limited, and 
therefore we have chosen not to include observational points 
at smaller radii 5 . To get a good measure of the scatter we gen- 
erated 10 4 mock observations, with normal distributed log- 
arithmic error bars of 0.22 dex on the total in each ring y,, 
giving an relative error of ~ 50% on y,- itself. For simplic- 

4 We have also tried to use Ym instead of Fiso. but it does not yield as good 
a fit. 

5 Including a central bin will for this cluster only decrease the error bars 
though. 




FIG. 9. — A set of mock observations based on the "Coma" cluster. The 
green diamonds, artificially scattered in the radial direction for visibility, are 
the mock observations, while the thin green lines are corresponding fitted 
profiles. The black line is the "Coma" SZ profile, with l-cr errors and the 
width of each ring overplotted. In the inset histogram is shown the distribu- 
tion of the reconstructed masses. The dotted line is the real mass, which is 
slightly offset compared to the average reconstructed mass. 



ity we have disregarded any correlation there may exist be- 
tween the different radial bins, and artificially have fixed the 
relative errors to a uniform level. The five rings are then fit- 
ted (see Fig. |9]l using Eq. ^ giving ro. The integrated SZ 
effect Ti8o for each mock observation is found by integrat- 
ing the fitted profile. Inserting ro and Y^o derived from each 
mock observation into Eq. (TBI) we find a reconstructed mass 
ofMi 8 o = 10.0 ± 1.4 x 10 14 M Q , or an 14% error on the recon- 
structed mass. There is a 4% systematic offset compared to 
the real mass of 10.4 x 1O 14 M . With five rings the error on 
the total signal is ~ \/5 50% = 22%. The mass goes roughly 
as T^gQ 7 , and we would expect roughly a 13% error on the 
mass, in agreement with what is found. 

4. DISCUSSION 

In this paper we have constructed a simple empirical model 
for the radial profile of the Sunyaev Zeldovich effect in groups 
and clusters of galaxies. The model has been motivated by, 
and validated against a mixture of simulations and observa- 
tions, and is characterised by only two parameters: The over- 
all normalisation, and a typical length scale related to the 
slope of the profile. It gives a very good fit to the simulated 
systems, and there is a tight relation between the parameters 
and the total mass of the system. Furthermore, the results are 
robust to the detailed gas physics employed in the simulations. 
This can be seen by considering the subset of the systems, the 
"Virgo" clusters, that are started from the same initial condi- 
tions, but simulated with different implementations of the gas 
physics (see Table 1). Because of the different gas physics 
there is an appreciable scatter in Y[&o, but still the mass (Mi go) 
of the clusters is well reconstructed (see Fig. [8j. 

The simulations are in good agreement with an observed 
average profile extracted from WMAP data for the outer part 
of the SZ profiles 0.2riso < r < 2r^o- Currently there are no 
published high resolution observations of the core SZ profile, 
but we can get an indirect measure by using temperature pro- 
files extracted from clusters observed in X-rays together with 
dark matter potentials from the simulations. Reconstructing 
the SZ profile, under the assumption of hydrostatic equilib- 
rium, we see that these "core profiles" are relative peaked, 
compared to the simulations. This is traced to differences in 



SZ profiles for clusters of galaxies 



7 



the temperature profile of the simulated systems, compared to 
what is observed using X-rays. We have to await future ob- 
servations of the resolved cluster cores, to determine to what 
extent this cusp in the central part is real, or a result of the hy- 
pothesis of hydrostatic equilibrium, which is known to be vi- 
olated in the centre of clusters of galaxies, where non-therma l 
processes such a s AGN heating dRovchowdhurv et al.ll2005l) . 

etal.l2l 



and cosmic rays dPfrommer et al.l2007l) can play an important 
role for the pressure balance. 

We stress that the model we have presented here is readily 
applicable to future observations. This will give a good 
proxy for the mass of the observed system. It will also help 
in reconstructing the full SZ profile from observations with 
low spatial resolution, to be used in conjunction with X-ray 



observations in the study of cluster dynamics. A fitting 
routine written in IDL for the profile that can be applied to 
observed or simulated data can be found at 

|http: 



77^ 



.phys . au . dk/~h augboe l/ software . shtml[ 
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